Thermal and electrical transport in single crystalline MgB2 
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We present the results of measurements of the thermal and electrical conductivity in the basal 
plane of single-crystalline MgB2 between 2 and 300 K. The analysis of the temperature dependence 
of the thermal conductivity gives supporting evidence for two different gaps on different sheets of 
the Fermi surface in the superconducting state. The zero-temperature values of the two gaps are 
Ai(0) = 5.3 ± 0.5 meV and A2(0) = 1.65 ± 0.2 meV. 

PACS numbers: 74.25.Fy, 74.70.-b, 74.25. Kc, 72.15.Eb, 



The recent discovery of superconductivity in MgB2 be- 
low an unexpectedly high critical temperature Tc of the 
order of 40 K initiated intensive studies of its physical 
properties [1] . Numerous results indicate that the super- 
conducting state of MgB2 is well described by the origi- 
nal BCS theory. Most experiments are compatible with 
a nodeless superconducting order parameter and imply 
that the electron-phonon interaction is responsible for 
the pairing of the electrons [2]. However, various types of 
experiments (see Ref. [3] and references therein), mainly 
using powder or polycrystalline samples, have given ev- 
idence for two gaps in the quasiparticle excitation spec- 
trum of this superconductor. Qualitative and especially 
quantitative experimental checks of this intriguing situa- 
tion on single crystalline material seem in order. 

Below we present results of measurements of the ther- 
mal conductivity k and the electrical resistivity p paral- 
lel to the basal a6-plane of the hexagonal crystal lattice 
of MgB2 as a function of temperature T between 2 and 
300 K. We demonstrate that our low-temperature results 
cannot be explained in terms of a single-gap function be- 
low the critical temperature, but are well approximated 
by assuming two different energy gaps in the supercon- 
ducting state. 

The investigated single crystal with dimensions of 
0.5 X 0.17 X 0.035 mm'^ was grown with a high-pressure 
cubic anvil technique as described elsewhere [4]. The 
resistivity was measured by employing a 4-contact con- 
figuration in the a6-plane. For the thermal conductivity 
measurements, the standard method with a constant uni- 
axial heat flow was used. The same contacts were used 
for measuring the voltage and the temperature difference, 
respectively. Additional measurements of p{T) and k{T) 
in constant magnetic fields if, oriented along the hexag- 
onal c-axis, were made as well. 

The electrical resistivity p{T) in zero magnetic field is 
presented in Fig. 1. The narrow (ATc = 0.1 K) super- 
conducting transition occurs at Tc = 38.1 K. In the tem- 
perature region between Tc and 130 K, p{T) may be well 
approximated hy p = p^+A'T^, where po and A' are con- 
stants. At temperatures below 50 K, p w po = 2.0 /if2 cm. 
A cubic p{T) dependence has often been observed in tran- 
sition metals and is associated with interband elcctron- 
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FIG. 1: In-plane electrical resistivity p{T) of hexagonal 
MgB2. The inset emphasizes the low temperature part for 
H = and 50 kOe along the c-axis. 



phonon scattering[5, 6]. The corresponding equation 

p{T)^ Po + AT^JsiSR/T), (1) 

where J„(z) — x^e^ (1 — e^)~^ dx, fits our data per- 
fectly. This is demonstrated by the solid line in Fig. 1. 
A similar fit to data from polycrystalline material is pre- 
sented in Ref. [7]. The fit parameters A and Qb. char- 
acterize the strength of electron-phonon interaction and 
the cut-off frequency Wmax of the phonon spectrum, re- 
spectively. The fit value of 8i? — 1226 K is in excellent 
agreement with the value of ?iwmax = 105 meV mea- 
sured in inelastic neutron scattering experiments [8, 9]. 
The magnitude of 8/j is considerably higher than the 
Debye temperature ~ 700 — 900 K, extracted from 
low-temperature specific heat measurements [10, 11, 12]. 
This discrepancy is not unexpected because electrical re- 
sistivity and specific heat involve different averages of the 
lattice frequencies [13]. 

The inset of Fig. 1 emphasizes p(T) close to the su- 
perconducting transitions in fields of iJ = and H = 
50 kOe. For = 50 kOe, the critical temperature, 
marked by the onset of the deviation from the constant 
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FIG. 2: Thermal conductivity vs temperature in the aft- 
plane of MgB2 in zero magnetic field and for H — 0.62 kOe 
and 50 kOe parallel to the c-axis. The solid line represents 
the typical zero-field behavior of k{T) measured on polycrys- 
talline samples (data from Ref [15]). 



resistivity in the normal state, is shifted to Tc = 17 K. 

The thermal conductivity data k{T) in zero magnetic 
field and for H = 0.62 and 50 kOe are presented in Fig. 2. 
The zero-field k{T) values are about an order of magni- 
tude higher than previously reported for polycrystalline 
samples [7, 14, 15, 16]. Also the temperature dependence 
of K is quite different from those earlier data. Instead of 
a monotonous increase with temperature we note a dis- 
tinct maximum of k(T) at T ~ 65 K. No anomaly in k{T) 
provides evidence for the transition at Tc. A pronounced 
slope change in k{T) is observed around T ~ 6 K, how- 
ever. 

The thermal conduction of a superconductor is usually 
provided by electrons (Kg) and phonons (Kph), such that 



+ Kph- 



(2) 



A separate identification of these two terms at arbitrary 
temperature is not straightforward. In the normal state, 
a convenient and often used way to estimate Kg is to em- 
ploy the Wiedemann-Pranz law, which relates the elec- 
trical resistivity and the electronic contribution to the 
thermal conductivity via 



Ke{T) = LoT/p{T), 



(3) 



where Lq = 2.45 x 10~^ W O is the Lorenz num- 
ber. The validity of this law requires an clastic scatter- 
ing of electrons and it is well established that Eq. (3) is 
applicable if the scattering of electrons by defects dom- 
inates. This is usually true at low temperatures, where 
p{T) ~ po- At higher temperatures but below Qd, 



Eq. (3) is usually invalid. In our case, the regime of appli- 
cability of Eq. (3) is limited to temperatures below about 
50 K (see Fig. 1). At temperatures > T > 50 K, Ke{T) 
provides about half of the total thermal conductivity. 

Upon decreasing the temperature below T^, the reduc- 
tion of the number of unpaired electrons leads to a de- 
crease of Kg and an increasing Kph. The overall behavior 
of k{T) in the superconducting state depends on the rel- 
ative magnitudes of Kg and Kph and also on the strength 
of electron-phonon interaction. Using existing theoreti- 
cal models, k{T) data may provide valuable information 
about the superconducting gap function. In the follow- 
ing, we restrict our analysis of k(T) to temperatures be- 
low 50 K, where the separation of Ke and Kph can reliably 
be accomplished. 

The electronic contribution to k{T) in the supercon- 
ducting state has been calculated by Bardeen, Rickayzen, 
and Tewordt [17]. In their model 
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as well as Fn{-y) = /o°°-2"(l + e''+y)-'^dz, y = 
A(T)//cbT, and A(T) representing the energy gap. Here 
Ke^n is the normal-state electronic thermal conductivity 
and in our analysis Ke,n = LoT/po. 

The phonon thermal conductivity was analyzed in 
terms of the Debye-type relaxation rate approximation. 
The corresponding equation is [18] 
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where w is the frequency of a phonon, t{lu, T) is the cor- 
responding relaxation time, v = 0£) (fcs/fi) (67r^n)~^/^ 
is the average sound velocity, n is the number density 
of atoms, and x = Tito/kBT. For our analysis, we use 
Bd = 750 K, as deduced from specific heat measure- 
ments [10, 11]. The total phonon relaxation rate 

r"^ = L/v + Blo* + CTlo^ exp(-ei5/6T) + Du 

+ELug{x,y) (7) 

can be represented as a sum of terms corresponding to in- 
dependent scattering mechanisms. The individual terms, 
dominating in different temperature intervals, introduce 
phonon scattering by sample boundaries, point defects, 
phonons, dislocations, and electrons, respectively. The 
constants L, B, C, b, D, and E are a measure for the 
intensity of corresponding phonon relaxation processes. 
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The function g{x,y) is given in Ref. [17]; we do not re- 
produce its rather complex form here. 

For H = 50 kOe the phonon contribution «;ph(r) in the 
normal state (above 17 K), extracted from k{T) and p{T) 
using Eqs. (2) and (3), is identical to Kph{T) calculated 
for if = above T^. We thus conclude that magnetic field 
effects on Kph are negligible. Based on the assumption 
that phonon scattering by both defects and phonons is 
the same in the normal and in the superconducting state, 
we estimate the temperatiuc region where the number of 
unpaired electrons is too small to produce significant con- 
tributions to both phonon scattering and heat transport. 
Taking into account the lowest previously claimed value 
of A(0) = 1.7 meV [19], the electronic quasiparticles are 
eff'ectively excluded from heat transport below 4 K. 

The following analysis procedure is based on the as- 
sumptions discussed above. In a first step, the values 
of the parameters L, B, C, b, D, and E, related to the 
processes of phonon relaxation in the normal state, were 
obtained by fitting the Kph{T) data in H = 50 kOe above 
17 K (for this region g{x,y) = 1) to Eqs. (6) and (7) as 
well as k{T) below 4 K where g{x, y) ~ and Kg = 0. 
These fits are shown in the inset of Fig. 3. The val- 
ues of the parameters, common for both fits, are L = 
2.6 X 10-4 m, B = 1.07 x lO"-'' C = 3.0 x lO^i* s K"i, 
h = 6.2, D = 2.1 X 10-^ and E = 7.5x 10"^ We refrain 
from a detailed discussion of the parameters related to 
defect, boimdary, and phonon- phonon scattering, and in 
the subsequent analysis we consider them as serving as a 
background. Our primary interest is the electron-phonon 
interaction represented by the parameter E. The value 
of E obviously fixes the shift between the two broken 
lines in the inset of Fig. 3. We note that the value of 
E is quite robust with respect to the choice of different 
phonon scattering terms in our fitting procedure. 

In the essential step of the analysis, we aimed at es- 
tablishing the temperature dependence of the gap func- 
tion A(T), using Eq. (7). The immediate result of this 
analysis is that a single gap function with the temper- 
ature dependence given by the BCS theory [20] is not 
adequate for MgB2. The significant reduction of the low- 
temperature thermal conductivity by a rather weak mag- 
netic field of 0.62 kOe (see Fig. 2) implies a strong in- 
teraction between phonons and the quasiparticles in the 
cores of the vortices induced by the magnetic field. In 
case of siich a strong clcctron-phonon interaction and 
the formation of a single BSC-type gap, a significant 
change of the slope of k{T) at is expected, as illus- 
trated in Fig. 3; where the dotted line represents the 
calculation of k{T) for A(0) = UQUbTc and A(r) as 
tabulated in Ref. [20]. Hence, for the majority of quasi- 
particles involved in phonon scattering below Tc, the 
gap opens much more slowly with decreasing tempera- 
ture than expected from the standard BCS prediction. 
Our fitting procedure is much more successful if we as- 
sume the formation of two energy gaps, as has previously 
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FIG. 3: In-plane thermal conductivity vs temperature in 
zero magnetic field. The solid line represents the best fit to 
the data assuming two gap functions (see text). The dot- 
ted curve results from a standard single-gap calculation with 
A(0) = 1.76A;b Tc. The inset represents the phonon thermal 
conductivity Wph in the regime of negligible phonon-electron 
scattering {H = 0, open circles), and the open squares repre- 
sent Kph(T) at temperatures where the electron-phonon scat- 
tering is that of the normal-state {H = 50 kOc, T > Tc). The 
solid lines represent the fits described in the text. 



been claimed in reports of other experimental studies of 
MgB2. Thus, we consider two subsystems of quasiparti- 
cles with gaps Ai and A2, different parameters Ei and 
E2 of phonon-electron scattering, and separate contribu- 
tions Kei and Ke2 to the heat transport. The normal-state 
parameter E is then given by the sum E = Ei + E2. In 
the superconducting state, the last term in Eq. (7) adopts 
the form 

Euj [(1 - a)g{x, yi) + ag{x, 2/2)] , (8) 

where yi = Ai{T)/kBT and y2 = A2{T)/kBT. The 
parameter a = E2/(Ei + E2) characterizes the relative 
weight of phonon scattering by quasiparticles in the two 
subsystems. Correspondingly, 

Ke = — [(l-<5)/(2/i)+(5/(y2)], (9) 
Po 

where the parameter 6 = k^i / and Kg = Kei + Ke2 ■ The 
fit is shown in Fig. 3 as the solid line. 

The resulting values of the fit parameters are Ai (0) = 
5.3 ± 0.5 meV, A2(0) = 1.65 ± 0.2 meV, a = 1 ± 0.02 
and S = 0.33 ± 0.02. The ratio 2Ai(0)/(fcBTc) = 3.22 
is close to the weak-coupling BCS value of 3.52. The 
second gap, however, is considerably smaller. Our iden- 
tification of two gaps is in qualitative and quantitative 
agreement with conclusions from results of specific heat. 



4 



point-contact spectroscopy, and photocmission experi- 
ments (see, e.g., Table I in Ref. [3]). The two gaps 
in MgB2 are, by most authors, associated with differ- 
ent sheets of the Fermi-surface. According to density- 
functional calculations of Liu et al. [21], the smaller and 
the larger gap are associated with 3-dimensional sheets 
and 2-dimcnsional tubes, respectively. However, the first- 
principle FLAPW band calculations of Hase and Yamaji 
[22] claim the opposite. Our result of a = 1 indicates 
that the dominant part of the interaction between quasi- 
particles and low-frequency phonons is provided by that 
part of the electronic excitation spectrum experiencing 
the smaller gap. The parameter E in Eq. (7) can be 
expressed as 



E : 



dof 
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where m* is an effective electron mass, E'def is the de- 
formation potential and Dq is the mass density. The 
effective masses of quasiparticlcs on different parts of the 
Fermi surface of MgB2 do not differ much [23]. It has 
been shown [24] that the interaction between the holes of 
the cylindrical a-bands and the in-plane bond-stretching 
optical mode E2g with a frequency of 470 cm~^ is pro- 
vided by an ultrastrong deformation potential. If E^ef 
is also strong for low-frequency acoustic phonons of the 
same polarisation, then our result implies that it is the 
holes of the 2D sheets of the Fermi surface that experi- 
ence the smaller energy gap, a coimterintuitive result but 
in agreement with the conclusions of Ref. [22] . Consider- 
ing the contributions of the two subgroups of quasipar- 
ticlcs to the heat transport, our result S = 0.33 suggests 
that the electrons experiencing the larger gap provide 
approximately 2/3 of total electronic heat conduction. 

In conclusion, we present a study of the thermal con- 
ductivity K in the a6- plane of a single crystal of MgB2. 
The analysis of the temperature dependence of k, pro- 
vides new convincing evidence for the formation of two 
superconducting gaps which are associated with different 
sheets of the Fermi surface. 
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